Open Charm and Beauty at Ultrarelativistic Heavy Ion Colliders 



Magdalena Djordjevic, 1 Miklos Gyulassy, 1 and Simon Wicks 1 

'Department of Physics, Columbia University, 538 West 120-th Street, New York, NY 10027 

(Dated: February 7, 2008) 

Important goals of RHIC and LHC experiments with ion beams include the creation and study of 
new forms of matter, such as the Quark Gluon Plasma. Heavy quark production and attenuation 
will provide unique tomographic probes of that matter. We predict the suppression pattern of open 
charm and beauty in Au + Au collisions at RHIC and LHC energies based on the DGLV formalism 
of radiative energy loss. A cancelation between effects due to the y/s energy dependence of the high 
Pt slope and heavy quark energy loss is predicted to lead to surprising similarity of heavy quark 
suppression at RHIC and LHC. 
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Introduction. RHIC and LHC experiments involving 
nuclear collisions are designed to create and explore new 
forms of matter, consisting of interacting quarks, anti- 
quarks and gluons. One primordial form of matter, called 
the Quark Gluon Plasma (QGP), is believed to have ex- 
isted only up to a microsecond after the "Big Bang" . If 
this QGP phase can be created in the laboratory, then a 
wide variety of probes and observables could be used to 
diagnose and map out its physical properties. 

The striking discoveries |l| at Relativistic Heavy Ion 
Collider (RHIC) of strong collective elliptic flow and light 
quark and gluon jet quenching, together with the deci- 
sive null control d + Au data, provide strong evidence 
that a strongly coupled Quark Gluon Plasma (sQGP), is 
created in central Au + Au collisions at V200 AGeV with 
gluon densities 10-100 times greater than nuclear matter 
densities 0. While there has been considerable conver- 
gence on the theoretical interpretation of RHIC data, 
the experimental exploration of the sQGP properties be- 
yond the discovery phase has barely begun .41]. Future 
measurements of rare probes such as direct photons, lep- 
tons, and heavy quarks will help to more fully map out 
the sQGP properties and dynamics. 

Heavy quarks provide important independent observ- 
ables that can probe the opacity and color field fluc- 
tuations in the sQGP produced in high energy nuclear 
collisions. In this letter, we present predictions of open 
charm and beauty quark suppression that can be tested 
at both RHIC and the future LHC facilities. Together 
with the already established light quark and gluon jet 
quenching and collective elliptic flow, a future observa- 
tion of a reduced heavy quark suppression (as compared 
to the observed pion suppression) could strengthen the 
current case for sQGP formation as well as test the evolv- 
ing theory of jet tomography 0. 

The prediction of D and B meson suppression pattern, 
in principle, requires theoretical control over the inter- 
play between many competing nuclear effects |(j that can 
modify the p± hadron spectra of heavy quarks. To study 
the high p± (p± > 6 GeV) heavy quark suppression, we 
concentrate on the interplay between two most important 
effects, i.e. jet quenching [£],0 and energy dependence of 



initial pQCD heavy quark p± distribution. In addition, 
we explore a range of initial conditions at LHC based on 
extrapolating RHIC data and based on Color Glass 
Condensate effective theory 0. We note that, for lower 
p± < 6 GeV spectra nonperturbative effects neglected 
here, for example collective hydrodynamic flow, quark 
coalescence and the strong gluon shadowing in the initial 
CGC state, may become important 0. 

Theoretical framework. To compute the heavy quark 
meson suppression we apply the DGLV generaliza- 
tion ^3 °f the GLV opacity expansion Q to heavy 
quarks. We take into account multi-gluon fluctuations as 
in ^l|. To apply this method, we need to know the fol- 
lowing: 1) Initial heavy quark p± distribution, 2) Differ- 
ence between medium and vacuum gluon radiation spec- 
trum and 3) Heavy quark fragmentation functions. 

The initial heavy quark p± distributions are computed 
using the MNR code [l^. As m tne El> we assume 
the charm mass to be M c = 1.2 GeV and beauty mass 
M& = 4.75 GeV. We assume the same factorization and 
renormalization scales as in |l3j |. For simplicity, we have 
concentrated only on bare quark distributions (< k\ >= 
GeV 2 ), and the runs were performed by using CTEQ5M 
parton distributions. 

Fig. ^ shows initial p± distributions for D and B 
mesons. By comparing p± distributions at RHIC and 
LHC case we see that RHIC distributions have signifi- 
cantly larger slope than the LHC ones. Since the sup- 
pression is sensitive to the slope of quark initial p± dis- 
tribution, the decrease in the p± slope with the increase 
of collision energy will have the tendency to lower the 
suppression from RHIC to LHC. 

Additionally, comparison between full and dashed 
curves on Fig^shows the variation of D and B meson p± 
distributions using two different types of fragmentation 
function. Full curves show meson spectrum obtained by 
using (^-function fragmentation, while dashed curves show 
meson spectrum obtained using the Peterson fragmenta- 
tion [Tj|- Though the choice of fragmentation function 
can lead to the order of magnitude difference in the ab- 
solute p±, we see that slopes of the curves remain quite 
similar. Therefore, we expect that the final suppression is 
insensitive to the choice of fragmentation functions. This 



2 




5 10 15 5 10 15 

Px [GeV] [GeV] 



FIG. 1: Initial p± distributions are shown for D (left figure) and B mesons (right figure). Lower (upper) curves 
correspond to RHIC (LHC) case. Solid curves are computed by assuming 5-function fragmentation while dashed curves 
assume Peterson fragmentation |l4j|. For D (B) mesons we used e = 0.06 (e = 0.006) 



conclusion is confirmed in Fig. 0] below, difference of less 
than 0.05 in the nuclear modification factor Raa is found. 
(Raa is the ratio of the observed yield in A + A divided 
by the binary collision scaled yield in p + p.) Therefore, 
for clarity, we show most results only for (5-function frag- 
mentation for both charm and beauty quarks. 

To compute the gluon radiation spectrum, we have to 
include (in general) three medium effects that control 
heavy quark energy loss. These effects are 1) The Ter- 
Mikayelian, or massive gluon effect [Hill, 2) Transition 
radiation pL7| which comes from the fact that medium 
has finite size and 3) Medium induced energy loss [lfl 
[T^ | , which corresponds to the additional gluon radiation 
induced by the interaction of the jet with the medium. 

In ^5 we will show that first two effects are not impor- 
tant for heavy quark suppression, since their contribution 
is less than 10% in the final result. Therefore, in this let- 
ter, we address only the medium induced gluon radiation 
spectrum which is given by fiof : 

dN™ = C F a s L [°° 2qVrfq 2 

dx " 7T \ g J (4|£)2 + ( q 2 +M2;r 2 +m 2)2 

f dk 2 6(2x(l - x)p x - |k|) 

X J ((|k|-|q|) 2 +M 2 ) 3/2 ((|k| + |q|) 2 +M 2 ) 3/2 

x { M 2 + (k 2 -q 2 ) ^~^~ TO U . (1) 
1' ; k 2 + M 2 x 2 +m 2 \ 

Here, k is the transverse momentum of the radiated 
gluon and q is the momentum transfer to the jet. M 
is heavy quark mass, p = 2(p/2) 1 / 3 is Debye mass, 

X g = 1 47 ^ a 2 p is mean free path ||, m g — p,/y/2 is gluon 

mass and E = \Jp\+ M 2 is initial heavy quark energy. 
We assume constant as = 0.3. For central collisions 
we take L = R x = R y = 6 fm, and assume that p 
is given by (1+1D Bjorken longitudinal expansion [Tflj , 
p = dN g / dyr-K L 2 , where is gluon rapidity density, 
and r is proper time. 



The energy loss was computed for both 1+1D Bjorken 
longitudinal expansion and using an effective average p 
approximation, where we replace r by < r >= -j. Since 
both procedures produce similar results, in this letter we 
present only on the computationally simpler (average p) 
results. 

We note that in Eq.(l) was set to k max — 2x{l — x)p± 
instead of k max — xE used in ■ Numerically, there is 
a 20% theoretical uncertainty in Raa due to the different 
reasonable choices of kinematical bounds. 

Heavy quark suppression at RHIC and LHC. In this 
section we compare suppression at RHIC and LHC as 
a function of momentum, collision energy and gluon ra- 
pidity density dependence. In Fig. we show Raa(j>±) 
for both charm and beauty quarks corresponding to D 
and B mesons in S fragmentation. For estimates of LHC 
initial conditions, we consider two cases: the PHOBOS 
extrapolation pfl (where gluon density is projected to be 
approximately 60% higher than at RHIC), and the CGC 
prediction p| (where the initial gluon density is predicted 
to be ~ 3 times higher than at RHIC). For charm quark 
we see that there is surprising similarity of Raa{p±) be- 
tween RHIC and LHC case, if PHOBOS extrapolation 
in gluon density is assumed. The similarity in suppres- 
sion between these results comes from the fact that, at 
LHC, the enhancement in energy loss (due to the larger 
gluon density), is mostly compensated by the decrease of 
the heavy quark distribution slopes. A slightly greater 
suppression is obtained with CGC estimate of the initial 
gluon density, which leads to larger energy loss. 

By comparing the charm and beauty suppressions on 
Fig. [21 we see that significantly less suppression is ex- 
pected for beauty than for charm quarks. This is due 
to the following two reasons 1) from Fig. ^ we see that 
beauty px distributions have significantly smaller slopes 
than the charm ones, and 2) due to dead cone effect |20| . 
the beauty energy loss is much smaller than charm en- 
ergy loss, as shown on Figs.l and 5 in [lOj ]. This explains 
in large part why no significant suppression was observed 
for pj_ > 2 GeV single electrons at RHIC [2l|. In this 
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FIG. 2: The suppression ratio Raa as a function of p±is 
shown for charm (lower curves) and beauty quarks (upper 
curves). Full curves correspond to RHIC case (yfs NN = 200 
GeV), while dashed and dot-dashed curves correspond to 
LHC case (\/s NN = 5.5 TeV). Dashed (dot-dashed) curves 
correspond to PHOBOS \J\ (CGC |8j) extrapolation in gluon 
rapidity density. 

kinematic range there is significant beauty contribution 
to the single electron yields and that component is essen- 
tially unquenched. Cronin and possibly collective flow 
effects in this low p± < 6 GeV region also may play a 
role. 

According to Fig. [2 we expect similar results for single 
electron suppression at both RHIC and LHC, i.e. we 
predict no significant suppression of single electrons at 
moderate pr at LHC as well. 

Our next goal is to study how the suppression is chang- 
ing as a function of collision energy. For that purpose we 
fix the p± of quark jet to 10 GeV and look at Raa{\/s) 
as shown on Fig. [3] We see that, if gluon density extrap- 
olates according to PHOBOS, than the RHIC « LHC 
conclusion from Fig. [2 is not a coincidence. It rather 
seams that, in this case, the high p± charm quark sup- 
pression is essentially independent on the collision en- 
ergy. In addition, the slight beauty suppression decreases 
as the collision energy increases. Therefore, we see that 
in PHOBOS extrapolation case, the 60% increase of the 
gluon density (and equivalently the increase in the en- 
ergy loss) is not enough to compensate the decrease in 
the p± slope. 

Slightly different situation occurs in the case of CGC 
extrapolation in gluon density. In this case, at LHC, we 
can expect 20% higher suppression for charm quarks and 
constant suppression for beauty quarks. 

Therefore, the main conclusion following from Figs. [21 
and |3 is that no significant difference between the RHIC 
and LHC heavy quark suppression is expected. This re- 
sult is surprising. To emphasize this point, we show in 
Fig. |3 dot-dashed curves showing a hypothetical case in 



which we assume that only energy loss changes with col- 
lision energy, while heavy quark initial p± distribution 
remains unchanged and fixed to 200 GeV case. From 
this curves we see that, at LHC, the energy loss would 
lead to additional 0.1 decrease in Raa for both PHOBOS 
and CGC case. 

If we compare the suppression for PHOBOS and CGC 
case on Fig. |3J we see that at 5.5 TeV (LHC) the dif- 
ference of 1000 in gluon rapidity density leads to only 
« 0.1 difference in Raa- Since the is still unknown 

at LHC, on Fig.Hwe show R AA (^) for a 10 GeV D and 
B mesons. We see that both D and B meson suppression 
falls slowly with the increase of the initial gluon rapidity 
density. 

Conclusions. In this letter we predicted the nuclear 
modification factor Raa(pt, Mq, \/s, f° r charm and 
beauty quark production in central Au 4- Au reactions 
with yfs = 200 - 5500 AGeV. We predict a rather weak 
\fs dependence in this range due to the compensation 
of the increasing energy loss in the more opaque sQGP 
and kinematic reduction of the px slope. Of course, it 
is still straightforward to deconvolute these competing 
effects to determine the growth of initial density with a/s 
and therefore differentiate between different predictions, 
such as CGC, of those initial conditions. 

By comparing our heavy quark predictions to the sup- 
pression patterns for the neutral pions in Ref. Q (light 
quark and gluon case), we expect a striking difference in 
the suppression pattern between light and heavy mesons. 
This is because the much more strongly quenched gluon 
jets component of light hadrons does not play a role in D 
and B production. The light hadron quenching pattern 
is therefore expected to have a stronger collision energy 
dependence @. 

We expect a moderate D meson suppression Raa ~ 
0.5 ± 0.1 for the ^ w 1000 ± 200 inferred from tt°. A 
similar suppression is expected at LHC for 1.5 — 3 times 
larger . Our high p± > 6 GeV predictions are robust 
within our approach, and significant experimental devi- 
ations would pose serious challenge to the pQCD based 
theory of radiative energy loss in sQGP matter. Future 
D meson data on 200 GeV d+Au and Au + Au and even- 
tually at LHC will thus enable critical consistency tests 
of the theory and the tomographic inferences drawn from 
the observed jet quenching patterns. 
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FIG. 3: The suppression ratio Raa as a function of y/(s) is shown for 10 GeV charm (lower curves) and beauty 
quarks (upper curves). Left (right) figure corresponds to the PHOBOS (CGC) extrapolation in gluon rapidity density. 
Upper x-axis show the gluon rapidity density that corresponds to the y/(s) for both PHOBOS and CGC scenario. Full 
curves represent the case where both energy loss and initial quark p± distribution change with -\/(s). Dot-dashed curves 

correspond to the case where only energy loss is changing with y/(s), while initial quark p±_ distribution is fixed at 200 
GeV. 
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FIG. 4: The suppression ratio Raa as a function of gluon 
density is shown for D (lower curves) and B mesons (up- 
per curves). Solid curves are computed by assuming 5- 
function fragmentation while dashed curves assume Peter- 
son fragmentation |l4|. For D (B) mesons we used e = 0.06 
(e = 0.006) 13]. 
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